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Density correlations unambiguously reveal the quantum nature of matter. Here, we study cor¬ 
relations between measurements of density in cold-atom clouds at different times at one position, 
and also at two separated positions. We take into account the effects of finite-size and -duration 
measurements made by light beams passing through the atom cloud. We specialise to the case of 
Bose gases in harmonic traps above critical temperature, for weakly-perturbative measurements. 

For overlapping measurement regions, shot-noise correlations revive after a trap oscillation period. 

For non-overlapping regions, bosonic correlations dominate at long times, and propagate at finite 
speeds. Finally, we give a realistic measurement protocol for performing such experiments. 
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I. INTRODUCTION 

Correlations between density measurements in quan¬ 
tum gases are related to the underlying correlation func¬ 
tions of the quantum field. In cold-atom experiments 
which detect individual atoms extracted from a cloud or 
a beam, the correlations appear in the arrival times and 
positions of the detected atoms [H-Q. Atoms are more 
(bosons) or less (fermions) likely to be coincident than 
classical statistics would suggest. 

By imaging a cloud after expansion and inspecting the 
fluctuations in the detected light, it is possible to infer 
correlations in momentum space[9j. This technique has 
been used to probe quantum many-body states such as 
the bosonic Mott insulator in a lattice and fermion 
pairing!^. Destructive in situ imaging has revealed spa¬ 
tial correlation information about one-dimensional qua¬ 
sicondensates (through both pair [T^ and triple [l^ corre¬ 
lations) and Fermi gases [Tj, [l^- More sophisticated op¬ 
tical detection schemes are used to measure correlations 
in spin polarisation[T^ [T^ rather than number density. 

Spatial density correlations can be inferred from fluc¬ 
tuations about the mean densityfl^ ■ Fluctuations are 
related to temperature and to thermodynamic suscep¬ 
tibilities like compressibility via fluctuation-dissipation 
theorems [ 3 . As susceptibilities tend to diverge close 
to phase transitions (e.g. Bose-Einstein condensation), 
so the measurable fluctuations become larger. However, 
fluctuation-dissipation theorems give information about 
the magnitude of fluctuations but not their length or 
time scales[l^. Naraschewski and Glauberfl^ showed 
how the underlying first-order correlation function re¬ 
lates straightforwardly to observable density-density[2l| 
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(second-order) correlations under the strong assumption 
that the system being measured is in the grand canonical 
ensemble at thermal equilibrium. In this paper we adapt 
their theory to include temporal correlations. 

Many of the measurements of correlation functions in 
atomic gases were done by detecting atoms which had al¬ 
ready been extracted from an expanding cloud or beam. 
Other experiments measured correlations between spa¬ 
tially separated regions of a cloud at one given instant 
or fluctuations from one experimental run to another. In 
the recent experiments of Guarrera et al. 1,0 atoms were 
extracted one by one from a single cloud and the tempo¬ 
ral correlations of their arrival times measured. There, 
the extraction of the first atom and its detection had no 
direct effect on the subsequent behaviour of the cloud: 
the correlations were present independently of the mea¬ 
surement, and the measurement itself was in some sense 
weakly perturbative [^. 

In contrast to previous experiments, we consider here 
density-density (not atom-atom) correlation measure¬ 
ments using light beams which pass through a single atom 
cloud at different times. The most general version of the 
measurement includes sampling the density at differing 
positions and times. These retarded correlations have 
not yet been observed in cold atoms. 

Correlations persist over thermal timescales like 
h/ksT at temperature T. A typical length scale is the 
thermal de Broglie wavelength: At = yhrWjMkrBT for 
atoms of mass M. Taking a cloud of ®^Rb at 200 nK 
the typical correlation times and lengths (without time- 
of-flight expansion) are 40 /iS and 0.4 ^m. Measurements 
can happen faster than 10 /iS [^ . with spatial resolution 
down to about 1 ^m, limited by photon collection and 
diffraction respectively. We therefore expect to be able 
to resolve correlations more easily in the time domain 
than in space. 

In this article we first state our method for calculating 
correlation functions for thermal bosons and then show 
how optical measurements of atomic density relate to the 
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correlation functions. We then specialise our discussion 
to the example of non-condensed Bose gases in harmonic 
traps. Finally, we discuss plausible experimental imple¬ 
mentations using weakly-destructive measurements. 


II. SECOND-ORDER CORRELATION 
MEASURES IN THE GRAND CANONICAL 
ENSEMBLE 


The purpose of defining this density correlation measure 
is to account for measurements over finite volumes and 
times, as will be shown in section Section PlBl 

Applying the definitions of G^^^and and the com¬ 
mutation relations, we note that ensemble average prod¬ 
uct of non-Hermitian operators can then be expressed in 
terms of the correlation functions: 






( 6 ) 


Following the formalism of Ref. , but retaining time 
dependence wherever possible, we start by defining the 
atomic-field annihilation operator for an atom at posi¬ 
tion r and time t as ^(r, t). Its Hermitian conjugate is 
the creation operator, and they follow the commutation 
relations for bosonic operators 




= 0 


( 1 ) 

( 2 ) 


Since G^^^ is real: 
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Re[Gg,^,^,,D(r,r',t-t') 


( 7 ) 


^(2) 

r' t t' 


For a non-interacting gas in a confining potential (with 
localised eigenfunctions) at thermal equilibrium, above 
the critical temperature for bosons, in the grand canoni¬ 
cal ensemble, it can be shown that(^ 


The function D depends on the Hamiltonian governing 
the system (see Appendix El . 

Density is an observable which is the expectation of 
the number density operator, n(r, t) = (h(r, t)) where 
n(r, t) = 'I'l(r, t)^(r, t). The quantum expectation value 
is denoted with angle brackets, and for readability we will 
use both notations n{r,t) and for expectation of the 
density at a specified position and time. 

We define normally-ordered, time-ordered (<' > t), 
first- and second-order atom-atom correlation func¬ 
tions [13, HH : 

= ( 3 ) 

cS'.tT = ( t')i>{v, t)) (4) 

The first-order correlation function (with equal time ar¬ 
guments) G*^^^ is the spatial representation of the one- 
particle density operator. Also, the density and correla¬ 
tion function are simply related: n(r,t) = G^^)(r, r, t, t). 
The correlation between pairs of atoms is found in G^). 

An obvious local, second-order, density correlation 
measure, as distinct to an atom-atom correlation func¬ 
tion, would be (hr.th-r'.t') — (h-r.t) (hr',t') . However, the 
operator in the first term is not Hermitian, and so this 
measure is not a good observable. We conjecture that 
the appropriate observable is simply the real part. 

We suppose that the measurements to be made are 
in some sense very weak, so that it does not matter in 
what order the measurements are made. The observ¬ 
able should then be symmetric under exchange of co¬ 
ordinates, {r,t} and {r',t'}. We note that when the 
variables are exchanged, the commutator in Eqn. @ is 
complex conjugated, and thus the real part is a correct, 
symmetric operator. Therefore, the correlation measure 
that we will use is: 

Gr,r',t,t' — 2 “t” ,t') (5) 


^( 2 ) _ ^( 1 ) ^( 1 ) 


-b 


G 


( 1 ) 
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where the second term is due to bosonic exchange symme¬ 
try. The correlation measure can be expressed in terms 
of the first-order correlation function G^^^ only: 


G. 




= Re 


G 


( 1 ) 


,D(r, r',t - t') -b 


G 


( 1 ) 


(9) 


The first term corresponds to counting (shot noise) num¬ 
ber fluctuations. The second term comes from quantum 
exchange correlations. 


A. The first-order correlation function for 
non-interacting, trapped, thermal bosons 


G(i) is conveniently expressed in a basis of the eigen¬ 
states of the trapping potential. First we write the field 
creation and annihilation operators in this basis set: 

'I'(r,<)=^u,(r)e-‘^^‘/'‘a, (10) 

all i 


where the state labelled i has energy e^, spatial vari¬ 
ation Ui{r) and annihilation operator fo. The sum 
runs over all states. The expectation {a\aj) = 6ijNi, 
with Ni being the thermal occupation number of 
the state, is given by the Bose-Einstein distribution, 
Ai = (l — . P = l/ksT is the inverse 

temperature, and ( = e^^ the fugacity with ^ being the 
chemical potential. It follows that the first-order corre¬ 
lation function is: 


G 


( 1 ) 


all i 


1 - C,e-P^- 


( 11 ) 
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For bosons, it turns out to be helpful to expand the 
thermal occupation factor in Eqn. m as a power se¬ 
ries of the fugacity, of the form x/{l — x) = 
Combining exponentiated terms, we find: 


GS.,,..=Ef 


all i 


The inner sum is in the same form as the representation 
of the propagator, or Green’s function[^. for the single¬ 
particle Hamiltonian H\ 

K{r,r\t,t') (13) 

all i 

= (r'|e-'^(*'-*)/''|r) 


Using this propagator with a complex time argument, the 
correlation function can be evaluated: 

OO 

G'S'.t,*' = E r', t, [t' - mp]) (14) 

fe=i 

At this point, we note that the expansion is effec¬ 
tively over increasing thermal occupation numbers, and 
that correlations with larger particle numbers will decay 
rapidly for non-condensed systems. 

The fugacity can be found by integrating the density 
over all space to give the total atom number: 


Ntotal = drrG 


( 1 ) 


(15) 


If the atom number is to be fixed in simulation, then the 
fugacity must be found so that this equation is satisfied. 


B. Correlation measurements over finite volumes 


spatial and temporal intensity profile /R.T(r,0 which is 
normalised such that f = 1 (dimensionless) 

at its spatial maximum [^. Photons arriving at detectors 
will be used to infer this atom number. 

A general second-order correlation measurement in¬ 
cludes a series of paired number measurements at two 
times, T and r' and with probe beams in modes centred 
at different positions R and R', as shown schematically in 
Fig.ffl This measurement corresponds to measuring the 
product of densities, which is the only kind of correlation 
that can be measured using near-resonant light [m. A 
good estimator for the correlations, the covariance, noted 
with the tilde to remind us that it’s spatially integrated 
and averaged over a series of realisations, is: 

GR,R',r.r' 

= JJ dtdt' JJ d^rd^r7R,T-(r,t)/R/_^/(r',t') (17) 


This finite-volume correlation estimator is the double in¬ 
tegral of the local correlation measure, Eqn. ®, weighted 
by the probe beam intensities. 



Imaging experiments measure the change of amplitude, 
phase or polarisation of a probe light beam due to the 
atom-light interaction[2l|. In the limit of small optical 
density, the change to the light beam is proportional to 
the average atomic density seen by the beam as it passes 
through the cloud, necessarily integrating along line of 
sight. There will also be a finite resolution transverse to 
the beam due to the optical geometry, ultimately limited 
by diffraction. Therefore, any optical measurement of 
density covers a finite volume. Real measurements also 
take finite times to count sufficient photons. This finite 
detection volume and time will tend to smear out and 
reduce measured correlations. 

The expectation of the time-averaged number of atoms 
in the probe beam during a measurement is [2^: 

E(A^R,r) = J dt J d^r {hr.t) lRrir,t) (16) 

with R a label for the central position of the probe beam 
for measurements around time r. The probe beam has a 


FIG. 1. A correlation measurement requires at least two atom 
number measurements at different times and places. Light 
beams pass through the atom cloud, and the number of atoms 
in the interaction volume shows up in the amplitude or phase 
properties of the light. 


This can in turn be expressed in terms of as in 
Eqn. ([9]). We find: 


Cn,w,r,T'= 11 dtdt' 11 d-^r d-^r' In,T(r,t)In' ,t')x 

2 ' 


Re 


G^'^Ct.t'D(r,r',t-t') 


r' t t' 


■ (18) 


It is common to work with correlation estimators which 
are normalised by the standard deviations of each of the 
component measurements. In this case, we choose to ap¬ 
proximate the standard deviations of measured numbers 
with Poissonian statistics in the limit of large number: 
Var (N) = E (N), where Var (A) is the sample variance. 
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The normalised correlation estimator is then: 


CR,RCT,r' = 


C'R,R',r,T' 

VE (A^r,,)E (iWO 


(19) 


III. LOCALISED CORRELATIONS IN BOSE 
GASES IN HARMONIC TRAPS 


We now turn to the specific case of thermal bosons 
in a harmonic trap. Our task is to evaluate the local 
correlation estimator, The Hamiltonian for an 

atom of mass M in a trap of frequencies and 

isH = -{h^/2M)V^ + + nlz^)/2. The 

eigenstates uimn(j) are the Hermite polynomials with en¬ 
ergies eimn = h(l^lx + mfly + nriz), with zero-point en¬ 
ergy subtracted for ease of notation. Since the Hamil¬ 
tonian is a sum of three terms, one for each dimension, 
both the eigenstates and the propagator are separable: 

= Kx{x, x', t, t')Ky{y, y', t, t')K^{z, z', t, t'). 

The propagator along the x direction, for the fcth term 
in the expansion of the correlation function, is given 

byH: 


{x, x',t, [t' — ihk/3]) = 

I mhI 


( 20 ) 


27ri?isin — t—ikhp)] 

iMVtx + x^) cosliixit' — t—ikhp)] — 2x'x^ 


exp 


2h sin[Ha; {f — t—ikhj5)] 


With a real argument this propagator oscillates, so there 
will be damped oscillations in the correlation function 
from the complex argument. 

The correlation function is then a sum over products 
of partial propagators: 

OO 

n (21) 

fc=l s=x,y,z 

Note that this does not mean that the correlation func¬ 
tion is separable. 

By numerical evaluation of the fugacity constraint 
Eqn. m- 


Ntotal=Y.^' n 

fc=l s=x,y,z 



it is possible to choose the atom number, and then to 
iteratively estimate the fugacity C to make this equation 
consistent. 


A. Two-position correlations at one time 

We now inspect the behaviour of the first-order corre¬ 
lation function for a cloud of 20 000 atoms of ®^Rb at 



FIG. 2. Un-normalised spatial correlation function for 
a cloud of 20 000 atoms of ®^Rb at 500 nK in a trap of 
1000 X 1000 X 20 Hz, along the centre line of the trap at the 
time zero. Since the phase-space density is 0.35, large parti¬ 
cle number occupancies (k) are not expected, and the series 
Eqn. (I21II from fe = 1 to A: = kmax converges rapidly. The 
fugacity which gives the correct particle number is ^ = 0.34. 
The typical correlation length is the de Broglie wavelength, 
which is around 265 nm in this case. 


500 nK in a cigar-shaped trap of 1 kHz x 1 kHz x 20 Hz, 
along the centre of the trap at time zero. In Fig. [21 we see 
that the series Eqn. (HID converges rapidly for tempera¬ 
tures significantly above critical. Each term in C,^ cor¬ 
responds to a thermal occupation number of the ground 
state by k particles. Since the temperature is above crit¬ 
ical for this example, the probability of thermal occupa¬ 
tion of large numbers of particles in the ground state is 
very small. The correlation function is significant over 
length scales around the de Broglie wavelength. 

For further discussion, we define the normalised corre¬ 
lation functions: 


- 


^(1) 


/^(l) ^(1) 

,r' ,t',t' 




^(2) 

'-^r r' t t' 


fy(k) ^( 1 ) 


(23) 

(24) 


which are independent of the number of particles in the 
trap and start at = 1 and = 2 for t = t' = 0 and 


B. Two-time correlations at one position 

Turning our attention to the normalised, two-time cor¬ 
relation function, Fig.|3]shows a decay of the two-particle 
correlations on the time scale Tc = hp (inset). The cor¬ 
relations revive and decay as particles oscillate in the 
confining potential, returning partially after the half the 
shortest oscillation period, i.e 500 /iS. The correlations 
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return completely on the longest time scale in the system, 
half the axial oscillation period (25 ms in this case). Any 
two particles will separate and then return to their orig¬ 
inal position after half an oscillation in the trap, so the 
two-particle correlation function revives. The real part 
of is negative for most of the time between the begin¬ 
ning and the first revival. This phase shift could in prin¬ 
ciple be observable in an two-path interferometer based 
on extracting atoms from the cloud, and then relatively 
delaying one of the paths. Beats between the different 
trap frequencies appear in the correlation function. 



FIG. 3. Normalised two-time correlation functions (real 
and imaginary parts) and — 1 for the same atom cloud 
as Fig. [2] at the centre of the trap. Inset: same figure but at 
shorter times. The correlation function decays on the thermal 
time scale Tc = fi/3 = 15 gs. Correlations revive after one half 
oscillation of atoms along the tighter axis, 500 gs= 32.7 x Tc 
here. 



t'! 


FIG. 4. Normalised two-time correlation function g^^^ for the 
same atom cloud as Fig. but offset from the trap centre 
along the short direction by a distance xq. We see that the 
revival goes from twice per oscillation period to just once as 
position is moved away from the trap centre. 


Inspecting at a fixed position but off-centre of the 
trap, Fig. 0] shows that the revival changes from occur¬ 
ring twice per oscillation period to just once. The dis¬ 
tance off-axis required for this effect is approximately the 
de Broglie wavelength. 


C. Two-time, two-position correlations 



FIG. 5. Normalised two-time correlation function g^^^ for the 
same atom cloud as Fig. [2] as a function of both time t' and 
position, with first point at the centre of the trap and time 
t = 0, the second off centre at x'. We see how correlations 
take a finite time to reach the second position. 

Fig. shows the temporal evolution of g^^^ when the 
measurement at the second time t' is done away from 
the trap centre. The correlation function peaks at later 
times at points further away from the centre. We can 
therefore define a speed at which the peak of the correla¬ 
tion function travels. It can be calculated by combining 
the typical correlation time scale Tc = hjS with the de 
Broglie wavelengt h At to make the correlation velocity 
Vc = \j2TrkBT/M. We see from figure Fig. [5] that corre¬ 
lations travel at about 2.5 x Vc- 


IV. MEASUREMENTS OVER FINITE 

VOLUMES IN OBLATE HARMONIC TRAPS 

We now apply the general form of Section IIIBI to the 
specific case of harmonically trapped Bose gases probed 
by Gaussian-shaped beams. The co-ordinate system is 
defined in Fig. [TJ We make the approximation that 
the Rayleigh range of the probe beam is much longer 
than the thickness of the cloud in the direction of prop¬ 
agation. This approximation makes it possible to sep¬ 
arate the detection response function, into components 
dR,r(r,<) = 

Measurements can be performed quickly compared 
to the typical correlation time scales, hf3, so we 
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will approximate the temporal response function of 
the detection with an infinitesimal response time: 

^R,r(r, t) = — t). The separated factors of the 

spatial detector response functions are Ix{x) = 1, 
Iviu) = and Iziz) = where 

R = {X, Y, Z) and wq is the waist of the probe beam. 

We consider a pancake shaped atomic cloud, which is 
thin along the line of sight in order to minimise the loss of 
coherence through integration along the probe axis; the 
other axes are kept loose. The trap contains 10® atoms 
at 100 nK and has trap frequencies of 413 x 20 x 9 Hz, 


much as in Ref. 

The full expression for the correlation estimator with a 
finite detection volume, Eqn. m, requires a numerically 
intractable 8-dimensional integral. Since the detector re¬ 
sponse contains a J-function we can execute the time inte¬ 
grals. The separability of both propagators and detector 
response can be used to express the 6-dimensional spa¬ 
tial integral as a product of two-dimensional integrals, 
reducing the numerical complexity. 

Making use of the result of Eqn. dSll) , we can now 
numerically evaluate Eqn. (HU as: 


=<^^C'^Re 

ffe=l Ls=x,y,z 


n JJ ‘isds'Is(s)Is'(s')K(^Hs,s',T, [r'- i?ik/3])K(°)(s',s,T',T) 


(25) 


oo oo 


, k—1 1—1 s—x,y,z 


EE ^ ds' Is{s)Is'is')K^'^'>*{s,s',T, [t' - ihfc/3])R:W(s,s',r, [r' - ihl/]]) 


The first term is due to quantisation of the matter field 
into atoms (shot noise). This term is dominant around 
the peaks of D{r, r', r' — r) (see Appendix IXll . It is max¬ 
imised at r' = T = 0, where it reduces to 

OO p 

C'rV=EC'= n / (26) 

k—1 s—x,y,z'^ 

Is{s)Is'{s)Ki’^Hs,s,0,[-ihkl3]). 

This equation gives an upper bound for the shot noise 
contributions. The second term in Eqn. (113 is due to 
bosonic fluctuations and will be labelled C®®. It decays 
more slowly than the first term when increasing t' — r. 

A. Two-position correlations at one time 

The correlation measure taken for two short (compared 
to h/3) measurements of the number of atoms in the probe 
beam, at nearly equal times, but two positions is shown in 
Fig. 13 The number of atoms in the beam estimated from 
equation Eqn. (113 is 37 atoms. The measures are taken 
10 /rs apart since the numerical integrals do not con¬ 
verge for exactly equal times. The bosonic-enhancement 
of fluctuations is up to 5 atoms. The bosonic correlations 
decay on the length scale of the measurement probe beam 
size, in this case 1.5 /xm. The bosonic-correlation length 
scale intrinsic to the atom cloud is At = 0.59 /xm here, 
much smaller than beam size. The shot-noise fluctuation 
at the origin is 24 atoms and therefore less than would 
be expected from Poissonian statistics, where variance 
equals the mean. The total mean atom number in the 
probe region is made of many atoms in the wings of the 
Gaussian probe beams, which are only partially detected, 
and hence fluctuations are reduced. 



FIG. 6. Un-normalised fluctuations C, for a cloud of 
10® atoms at 100 nK in a pancake-shape trap of 413 x 20 x 9 Hz 
as in Ref. [^, for two measurements at nearly equal times 
(just 10 fj,s apart). The phase-space density is 0.82. The scale 
over which correlations decay is the size of the probe beam 
used for detection, which is 1.5 pm in this case. For compar¬ 
ison, the de Broglie wavelength is 0.59 pm. The maximum 
shot noise fluctuation calculated with Eqn. dMI is 24 atoms. 
The mean atom number in the probe beam is 37, indicating 
that the fluctuations are effectively suppressed by the pres¬ 
ence of atoms in the Gaussian wings of the detection volume. 


B. Correlations at two times 

We now turn to normalised correlation estimators at 
two times as defined in Eqn. (ITOl) . Fig. [7] shows how the 
one-particle autocorrelation [first part of Eqn. (I7H1) ] de¬ 
cays faster than the bosonic correlations. Furthermore, 
the shot noise demonstrates anticorrelations as particles 
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which start in the middle move away. When they come 
back to the centre, a revival can be seen in the corre¬ 
lation function, roughly one full oscillation period along 
the tightest trap axis. 



FIG. 7. Normalised correlation estimator for the same cloud 
as Fig. [S] for two measurements at one position. The time 
scale over which correlations decay is the time a typical ex¬ 
citation takes to cross the probe beam used for detection, 
oc 1/%/T, and not the typical thermal timescale hfi = 76 /rs. 
At short times, shot noise (single-particle autocorrelation) is 
the dominant contribution to the correlation measure. Shot 
noise also exhibits a revival at about one oscillation period of 
the tightest axis of the trap. 

We also note that fluctuations become weaker as the 
temperature becomes much higher than the critical tem¬ 
perature for Bose-Einstein condensation. Furthermore, 
the time scale on which the correlations decrease is dic¬ 
tated by the time it takes atoms with the mean thermal 
speed to cross the detection region (again 1.5 p,m here). 
Since this speed increases with temperature, so the cor¬ 
relations die away faster at higher temperatures, with a 
typical timescale proportional to l/v^ and not the ther¬ 
mal coherence time scale, oc 1/T. 

When the two measurements at different times are also 
done at different positions, we see again that correlations 
take a finite time to move from the first point of mea¬ 
surement to the second. As shown in Fig. |51 when the 
probe beams are well separated {Z' > wq), the maximum 
of correlation is found at finite (not zero) time. This is 
an observable manifestation of the propagation of cor¬ 
relations at finite speed. For well-separated beams, the 
revival in correlations is not detectable. 


V. REALISTIC MEASUREMENT PROTOCOLS 
AND UNCERTAINTIES 

The determination of atom number will come with 
some unavoidable, fundamental uncertainties. We will 
now estimate how many measurements it will take to 



t' I flS 

FIG. 8. Normalised correlation estimator for the same cloud 
as Fig. [6] for two measurements at different positions as a 
function of the delay between measurements, t' . For the 
furthest separated measurement positions, the bosonic cor¬ 
relations propagate at finite speed. Negative autocorrelations 
cause the initial dip, and revivals are not observed. For com¬ 
parison, hp = 76 fis and l/fli, = 2420 fia. 


distinguish the bosonic part of the correlations from the 
shot noise. 

The uncertainty in the correlation measure, Var (C), 


assuming approximately Gaussian statistics, is[ 



1 




-Var(iVR,,)Var(fVR<^,0 


(27) 


where Var (Ar^^.) is the uncertainty in the measured 
number of atoms in the probe beam at position R and 
time r, and Nmns is the number of experimental runs. 

Optical detection of atoms always comes with incoher¬ 
ent photon scattering [ 2 ^, no matter what property of 
the light beam is measured. The more photons sponta¬ 
neously emitted, the more precise the measurement can 
be. However, spontaneous emission destroys the second- 
order coherence that we want to measure. We therefore 
have to limit the precision of each measurement, in favour 
of more, independent, experimental runs to reduce the 
statistical uncertainties. The measurement uncertainty 
per spontaneously emitted photon can be decreased by 
using small probe beams, so it pays to use high numeri¬ 
cal aperture optics. In a shot-noise limited measurement 
with a Gaussian beam of size wq the uncertainty in the 
inferred atom number is [ 2 ^: 


Var (Ar^,) = (28) 

ZcTqTIsc 

where n^c is the number of spontaneous emission events 
per atom and co is the resonant atom-light scattering 
cross section, which can be up to 3A^/27r for A the reso¬ 
nant optical wavelength. Note that this formula is only 
valid if the transverse size of the probe beam is signifi¬ 
cantly smaller than the atom cloud. 








































To detect the bosonic enhancement of the correlations, 
the variance on the correlations measure has to be less 
than the bosonic part of the correlation measure, the 
second term in Eqn. (1251) : 

Var (c) < C''" (29) 

Using Eqn. (HZl) and Eqn. (E51) we get an expression for 
the number of measurements required to measure the 
bosonic part of the correlation measure with a signal- 
to-noise ratio of one: 

(30) 

Taking the example of Ref. @ for experimental param¬ 
eters, using a 1.5 fim Gaussian-shape probe laser beam 
and limiting the number of photons scattered to 0.2 per 
atom, we find that a signal-to-noise ratio of 1 is achieved 
for about 750 experimental runs, for measurements at 
nearly equal positions and times. The cloud extends over 
a few times the size of the detection region, so it is pos¬ 
sible to measure correlations for several pairs of probe 
beams in parallel, pairs being separated by much more 
than intra-pair spacing. Several measurements can be 
made on one cloud, although atom number will reduce 
and temperature will rise. 

Only certain kinds of optical measurements are suit¬ 
able. The optical densities (exponents in the Beer- 
Lambert law) of the Bose gas clouds discussed are typ¬ 
ically larger than 1, so the detector response is not lin¬ 
ear in atom number and Eqn. HU) is not valid. Phase 
shift measurements (far off resonance) have a linear re¬ 
sponse, but will induce a net, localised, mean phase shift 
on the atoms due to the AC Stark effect. Such a phase 
shift will lead to mechanical effects on the atom cloud, 
such as heating or sound waves, which might disguise 
the atom-atom correlation signals. We have already de¬ 
veloped a suitable two-frequency measurement technique 
which cancels this mean phase shift . 


VI. CONCLUSIONS 


gases, and for other kinds of potentials, e.g. box poten¬ 
tials or optical lattices. Correlation measures for these 
systems would provide unambiguous evidence of exotic 
phase transitions. 

One can expect that strong, projective measure¬ 
ments will take the ensemble out of thermal equilib¬ 
rium. Extending the theory beyond thermal equilibrium 
is possible using, for example, stochastic time-evolution 
techniques [131 • We would be very interested to read 
about the results of any such calculations. 
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Appendix A: Properties of the commutation relation 
at unequal times 

In Eqn. m we decompose the field operators into a 
basis set of the eigenstates of the Hamiltonian that gov¬ 
erns the system. The commutation relations for creation 
and annihilation of particles within specific states are: 

a*, a] 

[ai,%]=0 (A2) 

where the Sij is the well-defined Kronecker delta. Explic¬ 
itly writing out the creation-annihilation operator for the 
fields we find: 

['k(r,t),'kt(U,t')] = ^ u.(r)u*(r')e+'(^A'-d)A[a„a]] 

all i,j 

all i 

=D{r,r'— t') (A3) 



We have shown that the results of temporal and spatio- 
temporal correlation measurements in non-interacting 
thermal Bose gases reflect the underlying atom-atom cor¬ 
relation functions. We have explained how to measure 
these quantities optically, for Bose gases in harmonic 
traps in thermal equilibrium. Correlations unsurprisingly 
take time to travel finite distances. The revivals in the 
correlation functions are in practice probably too weak 
to be visible in finite-volume correlation measurements, 
especially in the presence of non-zero atom-atom scatter¬ 
ing [lii. 

The two-time correlation functions can, in principle, 
be calculated in other quantum gas systems such as 
Bose-Einstein condensates, interacting Bose gases, Eermi 


With this formula we can calculate the unequal-time 
commutation relations for any bound (discrete) system, 
and the formula can be generalised to unbound (contin¬ 
uum) cases. By inspection, the commutator is directly 
related to the propagator of Eqn. Q. 

D{r,r',t — t') = K{r',r,t',t) (A4) 

At equal times Eqn. (IMl) reduces to the completeness 
relation for eigenfunctions, which is guaranteed by the 
fact that the Hamiltonian operator is Hermitian. 

At unequal times, the result depends on the Hamilto¬ 
nian. For our case of a 3D harmonic oscillator, which is 
trivially separable, there exists an analytic formula for 
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the propagator, Eqn. (nnn- This shows that the com¬ 
mutator is non-zero even for unequal positions. A phys¬ 
ical picture, for bound states, is that an initial delta- 
function wavepacket spreads out over time. Starting with 
two such wavepackets, with different positions at differ¬ 
ent times, there will be some overlap between them due 
to this spreading. Hence, their field operators do not 
commute, except for special cases like revivals in a har¬ 
monic oscillator. In those cases, the commutator tends 
to a Dirac delta function function of position in the limit 
where either the time arguments are equal or separated 
a multiple of the characteristic oscillation period. 

In general, a measurement done at time t will 
modify the state of the system and therefore change 
D{r, r', t — t'). For the purpose of this article, however, 
we made the strong assumption that the measurement is 


only weakly perturbative and can therefore be neglected. 

We also note that the numerical implementation of 
Eqn. (I25|) has two particular difficulties. First, the in¬ 
tegration of complex numbers is not well handled by 
some low-level integration routines. Secondly, and more 
severely, the propagators tend to Dirac delta functions 
of position in the limit where the time arguments are 
equal (or separated by an integer number of harmonic 
oscillator periods). That poses no problem for the sec¬ 
ond term due to the thermal, imaginary-time component, 
but means that numerical integrals of the first part can 
be non-convergent, or even worse, biased. In Fig. [51 we 
have not presented numerical results close to the diver¬ 
gences of the commutator for that reason. We note that 
evaluation of Fig. [8] took about 80 hours of processor 
time. 
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